Asymptotic entanglement of two atoms in squeezed light field 
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The dynamics of entanglement between two - level atoms interacting with a common squeezed reservoir is 
investigated. It is shown that for spatially separated atoms there is a unique asymptotic state depending on the 
distance between the atoms and the atom - photons detuning. In the regime of strong correlations there is a 
one - parameter family of asymptotic steady - states depending on initial conditions. In contrast to the thermal 
reservoir both types of asymptotic states can be entangled. We calculate the amount of entanglement in the 
system in terms of concurrence. 
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I. INTRODUCTION 

Dynamical creation of entanglement by the indirect interac- 
tion between otherwise decoupled systems has been recently 
studied by many researchers mainly in the case of two - level 
atoms interacting with the common vacuum. The idea that dis- 
sipation can create rather then destroy entanglement was put 
forward in several publications lUl-H]. In particular, the effect 
of spontaneous emission on the destruction and production of 
entanglement was discussed |5-8]. When the two atoms are 
separated by a small distance compared to the radiation wave- 
length, there is a substantial probability that a photon emitted 
by one atom will be absorbed by the other, and the resulting 
process of photon exchange produces correlations between the 
atoms. Such coiTelations may cause that initially separable 
states become entangled. 

The case of two atoms immersed in a common thermal 
reservoir was also investigated ll9l-[l2l. As was shown in lfT2l . 
similarly to the vacuum case, the collective properties of the 
atomic system can alter the decay process compared to the sin- 
gle atom. There are states with enhanced emission rates and 
such that the emission rate is reduced. The important example 
of the latter is the antisymmetric superposition \a) constructed 
from energy levels of considered atoms. When the atoms are 
close to each other, this state is decoupled from the environ- 
ment and therefore is stable. In that case the asymptotic states 
of the system are parametrized by the fidelity F of the ini- 
tial state with respect to the state |a) and the temperature T 
of the photon reservoir. Moreover, the asymptotic states can 
be identified with thermal generalization of Werner states i.e. 
mixtures of the state \a) and Gibbs equilibrium state at the 
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temperature T. 

In the present paper, we consider the atoms interacting 
with photon reservoir in a squeezed state ifisll . In practice, 
squeezed light sources produce photon fields in multimode 
squeezed states, but here we assume broadband approxima- 
tion in which the parameters characterizing the photon field 
are constant over a sufficiently broad frequency range. The 
dynamics of atoms interacting with squeezed light was stud- 
ied by many authors (see e.g. review paper [ 14] and references 
therein). In the context of our studies, we mention the result of 
Palma and Knight 1 15] showing the existence of highly corre- 
lated asymptotic state and the analysis of cooperative behavior 
of atoms in broadband squeezed light in Ref. 1 16]. 

In this paper, we study the asymptotic entanglement of the 
system of atoms evolving according the master equation con- 
sidered by Tanas and Ficek [18] but we allow non-zero detun- 
ing between the atomic transition frequency and the caiTier 
frequency of the photon field. In the case of spatially sep- 
arated atoms studied in details in Ref. ifTsll . there exists a 
unique asymptotic state, but in contrast to the vacuum or ther- 
mal reservoirs, this state can be entangled. However, the pro- 
duced entanglement is maximal only when the atoms are in 
resonance with the squeezed photon field. Non - zero detun- 
ing significantly diminishes this production. When the atoms 
are close to each other, the dynamics of the system radically 
changes. As in the vacuum or thermal case, the antisymmetiic 
state again decouples from the reservoir, therefore is stable. 
The asymptotic states Pas depend on the initial fidelity F and 
parameters describing the reservoir, but non - zero detuning 
also modifies the matrix elements of Pas- We show that the 
asymptotic states can be expressed as a mixture of a separable 
Gibbs state and two pure entangled states: an antisymmet- 
ric state I a) and some symmetric superposition of ground and 
excited levels of the atoms. This realization of the asymp- 
totic state simplifies for zero detuning and minimum - uncer- 
tainty squeezed reservoir to the mixture of \a) and two - atom 
squeezed state ifTsll . Thus in that case, there are two linearly 
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independent stable pure states so that decoherence - free sub- 
space is two - dimensional 1 17]. 

Depending on the initial fidelity, some of the asymptotic 
states are entangled. We calculate the amount of the asymp- 
totic entanglement using the concurrence as its measure. We 
also show that for initial fidelity greater than some threshold 
value (depending on the properties of the reservoir and detun- 
ing), the asymptotic concurrence is non - zero. This property 
is analogous to the thermal reservoir case. But when the reser- 
voir is in a squeezed state, somehow unexpected result occurs: 
initial states with small or even zero fidelity become asymp- 
totically entangled. The possibility of production of entangle- 
ment starting from separable states with zero fidelity is very 
interesting. In this case the correlations present in a squeezed 
reservoir are transferred to the atomic system, entangling for 
example two atoms both in the ground state. But as before, 
large detuning between atoms and photon field, destroys this 
possibility. 



II. MODEL DYNAMICS 

Consider two-level atoms A and B with ground states |0); 
and excited states 1 1) ^ (j =A,B), interacting with the radiation 
field in a broadband squeezed vacuum state with the carrier 
frequency 03^. The parameters and M characterizing the 
squeezing satisfy 



M = \M\ e'^ and \M\ < y^N{N+l), 

where the equality holds for a minimum-uncertainty squeezed 
state. In the Markov approximation the influence of the reser- 
voir on the system of atoms can be described by the dynamical 
semi-group with the Lindblad generator llisll 



where 



L = -i[H,-]+Lu, 



E E ^jkoia^, (II.1) 

j=A,B i.k=A.B 



and 



^ j.k=A.B 

j,k=A,B 



j.k=A,B 

\ E Y;iA?(2oipo* -oiaip-po^ai 



■ i,k=A,B 



(11.2) 



Here 



In the Hamiltonian ( III, lb . cOo is the frequency of the transi- 
tion \0) j \ l)j (j — A, B) and = Q.^^ = ^ describes in- 
teratomic coupling by the dipole-dipole interaction. On the 
other hand, dissipative dynamics is given by the generator 
( III.2I 1 with parameters Jab satisfying 



JAA = JBB = To, Jab = Jba = T 



(11.3) 



In the above equalities, yo is the single atom spontaneous 
emission rate, and y = G{rAB) Jo is the collective damping con- 
stant. In the model considered, G{rAB) is the function of the 
interatomic distance F^b, and G(?ab) is small for large sepa- 
ration of atoms. On the other hand, G{rAB) 1 when Vab is 
small (for more details see e.g. 1 19]). 

The time evolution of the system of atoms is given by the 
master equation 



dp 
dt 



Lp, 



(11.4) 



In the frame rotating with frequency (Oj, the master equation 
(III.4b becomes an equation with time indendent coefficients, 
and it may be written as 



dpi 
dt 



Lpi, 



(11.5) 



where 



with 



L=-i[H,-]+LD, 



H=^ oi+ ^jko'+ai, 5o = coo-co,, (II.6) 
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j=A.B j.k=A.B 



and 



LdPi^:: E Yi<r(l +A^) (2oip/oi - o^aip/ - p/O^oi 



■ j.k=A.B 

I- ^ E '^Jk^ (2o;p/oi - oio|p/ - p/oio 

j,k=A,B 



4 E Y;i^ (2o>/0^-oX+P/-P/<'+< 

j.k=A,B 

\ E Y;iA?(2oip/ol-oioip/-p/0*ai 



■ j,k=A,B 



(11.7) 



Notice that in the Hamiltonian ( III.6I ), detuning 5o can be arbi- 
trary. Only when the atoms are in resonance with the carrier 
frequency of the squeezed vacuum, 5o = 0. 

From now on we omit the subscript /. The master equation 
(111.5b can be used to obtain the equations for matrix elements 
of a state p of the system of two-level atoms with respect to 
some basis. To simplify the calculations one can work in the 
basis of collective states in the Hilbert space ® lfl9ll . 
given by product vectors 

k) = |i)A®|i>B, I^) = |o>a®|o)b, (11.8) 
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symmetric superposition 



and finally 



-^{\0)a®\1)b + \1)a®\0)b), (11.9) 



and antisymmetric superposition 



-j={\1)a®\0)b-\0)a®\1)b). 



(11.10) 



In the basis of collective states, two-atom system can be 
treated as a single four-level system with ground state \g), ex- 
cited state |e) and two intermediate states \s) and \a). From 
(IlI.SI l it follows that the matrix elements of the state p with re- 
spect to the basis |e), \s), \a), \g) satisfy the equations which 
can be grouped into decoupled systems of differential equa- 
tions. So the diagonal matrix elements and peg, satisfy 



dpe 



dt 

dpss 
dt 

dPaa 

dt 
dt 

dPeg 

dt 



(yo -y)A^P«« + iyo+y)Npss ~ 2yaNpee 

-l{Mpge+Mpeg), 

-(Yo+y)[(l +2A?)p,,- (1 +n)p,e-Npgg 

-Mpge-Mpeg], 

-{lO-y)[{l+2N)Paa-{l+N)Pe,-N Pgg 
+ Mpge+Mpeg], 

(yo - y) ( 1 + N) paa + (yo + y) ( 1 + A^) p.. - 2yo p,,^ 

-l{MPge+Mpeg), 

- (yo - y)p«« + (yo + y) m p,, - yM p^^^. 

-(yo(l+2A?) + 2/5o)p.g. 

(11.11) 



On the other hand, the elements Pae, Pag, Pse and Psg are con- 
nected by the following equations 



dpa 



dt 



dp 



ag 



dt 



dt 



dPsg 
dt 



y ( A' + ^ j - yo [2A^ + ^ ) + / (5o + i2) 



- (yo - y) Pga + (yo - y)Mp^„ - jMpag, 

y(N+^]-yo(2N+^)-i (8o - H) 



+ (yo - y) M pj,„ - (yo - y) ( 1 + A^) Pea - jM p„ 



y(A^+^)+yo(2^+^ ) -/(5o-i2) 



- (yo +y)Mp,, + {yo+y)Npgs - yMp,^, 
y(N+^-)+yo(2N+^-\+i (5o + il) 



+ (yo + y) ( 1 + A^) pes + (yo + y) m p^,. - yM p^^. , 

(11.12) 



dPas 

dt 



[yo(l+2A?)-2/i2]p„, 



(11.13) 



The equations for the remaining matrix elements can be ob- 
tained by using hermiticity of p. 

From the equations (111, lib it follows that similarly as in the 
case of reservoir in the vacuum state (see e.g. fl^) and ther- 
mal state |fT2ll . the system of atoms in the symmetric state \s) 
decays with the enhanced rate yo + y, whereas antisymmet- 
ric initial state \a) leads to the reduced rate yo — y. When the 
atoms are so close to each other that we can ignore the effects 
of their different spatial positions, we can put y = yo- In this 
limiting case of strongly correlated atoms (Dicke model) the 
state \a) is completely decoupled from the reservoir It can 
be also checked that the master equation (111.5b describes two 
types of time evolution of the system of atoms, depending on 
the relation between y and yo. When y < yo, there is a unique 
asymptotic state. This state was found in Ref. fl8^ for the spe- 
cial case of zero detuning. In general case we compute it in 
the next section. On the other hand, in the Dicke model case 
when y = yo, we show that there is a one - parameter family 
of nontrivial asymptotic states depending on the initial states. 



III. ASYMPTOTIC STATES 

A. Spatially separated atoms 

We start with the case of spatially separated atoms, when 
y < yo. Direct calculations show that in that case, there exists 
a unique stationary asymptotic state Pu, which in the canonical 
basis 

|1)a®|0)b, |0)a®|1)b, |0)a®|0)b 
has non-vanishing matrix elements 



(III.1) 



where for 



we have 



(111.2) 



ao 






_ CO 






P22 = 


P33 






mo' 






mo' 




bo 




zo 




do 




Pl4 = 




P44 = 




mo' 




mo' 




MO 



s §0 ^ y 
5 = —, y=— , 
yo yo 



Mo-(l+2A?)2 [(l+2A?)2+452] 
+ 4|M|2(f -(1+2A?)2), 



and 



flo=A^^ [(l+2A^)2-4|M|2 + 452] +|M|2f , 
co=NiN+l) [(l+2A?)2-4|M|2 + 452] +|M|2f , (III.3) 
do = {l+N)^ [{l+2N)^-4\M\^ + 45^]+\M\^f. 
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Moreover, 



bo = -2y|M|2, 

zo = -il+2N-2id)yM. 



(III.4) 



The state jllLlI ) in contrast to the analogous asymptotic state 
in the thermal reservoir, can be entangled and, as we show 
later, its entanglement crucially depends on the value of the 
normalized detuning 5 and the normalized damping constant 



B. Strongly correlated atoms 



When Y = yo, equations ( III. 1 11 ) - ( III.13I ) simplify and one 
can check that the solutions of ( |II.12t and ( III.13I ) asymptoti- 
cally vanish and the only contribution to the asymptotic states 
pas comes from Pee, Paa, Pss, Pgg and Peg. Note that in this 
case 



dPa 



dt 



where 



= 0, so p„„(f) ==p„„(0) 



F = (fl|p|a) 



is the fidelity of the initial state p with respect to the antisym- 
metric state I a). Hence the fidelity the asymptotic state Pas 
also equals to F and one finds that in the canonical basis the 
matrix of Pas has the same "X" form as in the case of the state 
(llll.lb . but with non - vanishing matrix elements given by 



pll = 


(1 


-F) 


a 
u 


P22 = 


(1 


~F) 


c 

2u 


P23 ^ 


(1 




c 

2u 


Pl4 = 


(1 


~F) 


z 
u 


P44 = 


(1 


^F) 


d 

1 



F 

2' 
F 

r 



(III.5) 



and P33 = P22- In the equations ( IIII.5I ) we have 

u ={l +2Nf {I +3N + 3N^ -3\M\^) 



-4(1+3A?+3A?2)5^ 



(III.6) 



and 



a^m^- [N{N+\)~\M\^] +|M|2 

c = (1 +2Nf [N{N+l) - |M|2] 

+ 2A'(A'+1)52, (III.7) 
d={l+2N) [l+N + 3{N{N+l)-\M\^)\ 

+ 2N [N{N+ 1) - |Mp] +4(1 +N)^5^ 
z = -{l+2N-2id)M. 



The asymptotic states Pas defined by (IIII.5I I exists for any ini- 
tial state, and for fixed parameters characterizing the squeez- 
ing depend on the initial fidelity and the normalized detuning 
8 = 5o/Yo of the electromagnetic field. When M = 0, we re- 
cover the case of standard thermal bath with playing the 
role of the mean photon number lfl2ll . 

To study the structure of the asymptotic states, we consider 
first the special case of a minimum-uncertainty squeezing and 
zero detuning of the radiation field. One can check that in that 
case, the matrix elements of Pas are given by 



Pll 


= {l-F] 


P22 


= (Pas)33 




F 


P23 


~"2' 


P44 


= {i-f: 


Pl4 


= ii-F: 



1+2N' 
_ F 
~2' 



l+N 
1+2N' 



(III.8) 



1+2A^ 



where 9 = + 7t. The asymptotic state given by dlll.SI l has a 
remarkable structure: it is a mixture 

p,s^{l-F)\N,Q)(N,Q\+F\a){a\ (111.9) 

of the pure state 



(III. 10) 

and the antisymmetric state \a). The state |A^,6) is known 
as two - atom squeezed state, and can be obtained from the 
ground state \g) — \0)a fX" |0)b by applying the atomic squeez- 
ing transformation S{t,), given by 



(III.ll) 



for the appropriate choice of the complex parameter t, lll5ll . 
This state is entangled and in the limit of maximal squeez- 
ing (A^ — > oo), it becomes a maximally entangled generalized 
Bell state. Notice also that \a) and |A^,9) span decoherence - 
free subspace for this specific system, as it was recently estab- 
lished in Ref. IitIi . 

In a general case the structure of Pas is much more involved. 
Define 



Frr — 



(III. 12) 



By a direct calculation we see that if F > Fcr, then 

pas = (l-p-9)pp + p|a)(fl|+?|¥)(¥|, (III-13) 

where 

.= (1 + -)^--, 
V u/ u 



u\J ad 
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The state Pp is a Gibbs state 



A. Entanglement of the asymptotic state Pu 



for the Hamiltonian Ha = Hq+ Hi with 



(III. 15) 



COi 



'3 ' 



j=A,B 

the inverse temperature 



R 1 1 ^ 

2(00 a 



and the frequency 



2(00 , c 
COl = : — — In ■ 



Ind/a y/^~\z\ 
Moreover, the pure state is given by 



(III. 16) 



(III. 17) 



l¥) 



10)^(8 |o)«+.'-%/^ 

\ a+d \ a+d 



\1)a<E>\1)b, (111.18) 



where (|) = arg z. 

The formula ( IIII.13I ) is a generalization of the equation 
( |III.9t as well as the corresponding representation of Pas by 
the thermal generalization of Werner states in the case of ther- 
mal reservoir lfT2l . Observe also that for F < Fa-, the asymp- 
totic state cannot be expressed as the mixture (|III.13l l but in 
contrast to the purely thermal case, the states Pas can be en- 
tangled even if F < F^r- We will study this problem in the next 
section. 



IV. ASYMPTOTIC ENTANGLEMENT 



Let us start with spatially separated atoms which have the 
unique asymptotic state Pu. Its concurrence is given by 



C(pu) = 2 max 



|zo|-co 1^0 1 - Vflo'^o 



MQ 



MO 



(IV.4) 



Analysis of this function in the general case of broadband 
squeezed reservoir is involved, so we focus on the case of min- 
imum - uncertainty squeezed states and consider (IIV.4t as the 
function of squeezed field intensity A^, for fixed values of pa- 
rameters y and 5. We plot this function in FIG. 1 for different 
values of detuning. It is evident that there is a range of values 
of mean photon number for which the asymptotic concur- 
rence is positive. Observe that the maximum of C(pu) appears 
for rather small values of A^ and the non - zero detuning dimin- 
ishes the production of entanglement. 
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FIG. 1: Entanglement of the state Pu as a function of N for Y= 0.85 
and different values of detuning: 6 = (dotted curve); 5 = 0.5 
(dashed curve) and 8=1 (solid curve). 



For the characterization of entanglement of the asymptotic 
state Pas we use Wootters concurrence li20ll defined for any 
two-qubit state p as 



C(p) = max ( 0, VX 



X4 , (IV. 1) 



where A-i > A,2 > A,3 > A-4 are the eigenvalues of the matrix pp 
with p given by 

p = a2(Xi02P 02(8)02, 

where p denotes complex conjugation of the matrix p. For the 
states in the "X" form, concurrence is given by the function 



with 



C(p) =max (0,Ci, C2), 

Cl =2(|pi4|-VP22P33), 
C2=2(|P23|-VP11P44)- 



(IV.2) 



(IV.3) 



B. Entanglement of the states Pas 

The properties of the concurrence of Pas as a function of 
initial fidelity can be studied in more details. Notice that for 
these states we have 



Define 



Fi — max 0, 



IfFi >0, then 



c-2\z\ 
c — 2\z\ — M 



Cl > for 0<F <Fi. 



On the other hand 

C2=2 



^ C F 



-{l-F) 



ad 



(IV.5) 



(IV.6) 



(IV.7) 
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Notice that if F < Fcr, then 



and 



Since 



^ C F 



„ / 2yad ~ c \ 2\/ad — c 
C2= — 1 . 



(IV.8) 



2\fad~ c 



1 <0, 



so 



C2 < when F < F„. 
Let now F > Fc^, then 

/ c + 2\/ad\ c + 2\/ad 
Ci = IH F . 



Define 



Fi 



c + 2^/ad 



(IV.9) 



(IV. 10) 



c + 2y/ad + u 
Form the equation ( |IV.9t we see that 

C2 > when F > Fi- 

On the other hand, direct calculations show that 

Fi > Fci and Fi < F2. 

Taking into account the above results we arrive at the con- 
clusion that depending on the initial fidelity F, the asymp- 
totic state Pas is entangled for all F e [0,Fi) U (F2, 1] (pro- 
vided Fi > 0) and separable for F G [Fi ,^2] (see FIG.2). The 
asymptotic concurrence reads 



C(Pas) 



Ci , 0<F <Fi 
C2, F2<F<\ 



(IV. 11) 



with Ci and C2 given by equations ( IIV.5I ) and ( |IV.9l l. respec- 
tively. This general result covers also the special cases of vac- 
uum reservoir where F2 — 0, thermal reservoir with F\~Q and 
F2 > and minimum-uncertainty squeezed reservoir, where 
Fi = F2. It is worth to stress that the creation of the asymptotic 
states with non-zero entanglement from the initial states with 
small or even zero fidelity is only possible when the reservoir 
is in a squeezed state. Let us discuss this point in more de- 
tails in a special case of atoms which are in resonance with 
minimum - uncertainty radiation field. In this case 



and 



Fi^F2 = 



C(Pas) ^ 



2^N{N+l) 



2y/NiN+l) + {l+2Ny 




with 



Co 



VN{N+1) 
1+2N 



(IV. 12) 



(IV. 13) 



(IV. 14) 



For all initial states with zero fidelity, we obtain pure entan- 
gled state dm. 101 ) with concurrence equal to Co. Notice that in 
the limit of maximal squeezing, this state becomes maximally 
entangled. For pure product states 



the fidelity is given by the formula 

F = \{i-\{<9m')., 



(IV 15) 



(IV. 16) 



so the zero fidelity corresponds, for example, to the case of 
two atoms prepared in the same initial states. This leads to 
the remarkable result: the interaction with squeezed reservoir 
will entangle two atoms which are initially in the ground state 
\g) — |0)a |0)b. The analogous phenomenon cannot occur 
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FIG. 2: Asymptotic entanglement versus fidelity for minimum - un- 
certainty squeezed reservoir with A' = 1 and detuning 5 = 0.8 
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FIG. 3: Asymptotic entanglement of initial states with F = as a 
function of N, for different values of detuning: 6 = (dotted curve); 
8 = 0.8 (dashed curve) and 8 = 2 (solid curve) 



7 



when the photon field is in the vacuum or thermal state. Notice 
also that for non - zero detuning, the asymptotic state is no 
longer pure and the production of stationary entanglement is 
less effective (FIG.3). 

V. CONCLUSIONS 

We have investigated the dynamics of two - level atoms in- 
teracting with the photon reservoir in a broadband squeezed 
vacuum state. Time evolution of the system crucially depend 
on the relative distance between the atoms. When the atoms 
are spatially separated, there is a unique asymptotic state, 
which can be entangled in contrast to the analogous asymp- 



totic state for the thermal reservoir In the case of small in- 
teratomic distance, there are nontrivial asymptotic states Pas 
which are parametrized by the fidelity F and the parameters 
N and M characterizing the squeezing. The states pas depend 
also on the detuning between atomic transition frequency and 
carrier frequency of the photon field. For the values of F 
above the threshold fideUty F2, the states pas are entangled. 
Non - zero entanglement can also occur for small values of 
F or even if F = 0, and this possibility is a unique feature of 
the squeezed reservoir. When the atoms are in resonance with 
the photon field and |M| = ■s/N{N+ 1), the asymptotic state 
corresponding to F = is the pure entangled state known as 
two - atom squeezed state. 
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